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1.0  INTRODUCTION 

Initial  reentry  and  space  vehicles  of  primary  interest  were,  in 
general,  axisymmetric.  A  capability  for  free-flight  testing  of  such 
configurations  having  force  and  moment  characteristics  linear  with  angle 
of  attack  hasj  existed  for  some  time,  both  in  wind  tunnels  and  in  aerobal- 
listic  ranges,  for  example.  Refs.  1  and  2.  With  the  imposed  restraint 
of  linear  force  and  moment  characteristics,  closed  form  solution  of  the 
differential  equations  of  motion  could  be  obtained  and  fitted  to  the  meas¬ 
ured  motion  histories  experienced  by  the  bodies.  In  such  a  fitting  proce¬ 
dure,  a  fit  of  the  measured  motion  history  is  obtained  by  adjusting  the 
unknown  parameters  in  the  motion  equations.  More  recently  the  develop¬ 
ment  of  numerical -integration  fitting  procedures  (see  Refs.  3  and  4),  in 
which  the  differential  equations  of  motion  are  used  directly  in  fitting  the 
measured  motion  histories,  has  permitted  handling  axisymmetric  config¬ 
urations  having  nonlinear  aerodynamic  coefficients.  However,  many  of 
the  current  and  proposed  vehicles  of  interest  have  appreciable  asym¬ 
metries.  Motion  histories  of  such  vehicles,  in  general,  cannot  be  de¬ 
scribed  adequately  with  the  standard  motion  equations  that  have  been 
used  satisfactorily  with  axisymmetric  bodies.  Hence,  asymmetric  vehi¬ 
cles  have  represented  an  area  in  which  the  testing  capability  was  inade¬ 
quate.  It  should  be  noted  that  asymmetries  as  used  here  refer  to  inertia 
asymmetries  and  to  primary  aerod3mamic  asymmetries  which  involve 
significant  differences  between  Cnq,  and  C^g  and  between  ICmal  and 

|Cjjg| .  This  is  in  contrast  to  publications  in  which  asymmetries  refer 
to  aerodynamic  trim  terms. 

This  report  presents  a  discussion  of  the  development  of  a  procedure 
(numerical-integration  fitting  approach  in  conjunction  with  revised  motion 
equations)  compatible  with  bodies  having  appreciable  aerodynamic  and 
inertia  asymmetries.  Further,  free-flight  experimental  data  for  elliptic 
cross-section  bodies  are  presented.  The  experiments  were  conducted 
in  the  von  K^man  Gas  Dynamics  Facility  (VKF),  AEDC,  1000-ft  Hyper- 
ballistic  Range  G. 

2.0  MOTION  EQUATIONS  FOR  ASYMMETRIC  BODIES 


The  conventional  motion  fitting  technique  as  used  for  many  years 
is  restricted  to  axisymmetric  bodies  and,  as  noted  previously,  is  re¬ 
stricted  to  bodies  having  linear  force  and  moment  characteristics  (see 
Ref.  5).  The  closed  form  solution  of  the  yawing-motion  differential 
equations  associated  with  that  technique  is  also  dependent  on  the  ratio  of 
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the  rolling  velocity  of  the  body  to  its  longitudinal  velocity  being  constant. 
Another  area  of  concern  in  that  fitting  technique  is  the  axis  systems  in¬ 
volved.  Although  the  experimental  yawing  motion  of  a  body  (measured  with 
a  range  shadowgraph  system)  is  defined  relative  to  a  fixed-plane  axis  sys¬ 
tem,  the  yawing-motion  equation  used  is  referenced  to  a  nonrolling  axis 
system  so  that  the  equation  can  be  written  in  closed  form.  A  sketch  of  the 
fixed-plane  axis  system  is  shown  in  Fig.  la.  In  this  system  0=0  and,  by 
definition,  the  y-axis  of  the  fixed-plane  axis  system  remains  in  the  (xo-yo) 
plane.  In  the  nonrolling  axis  system,  the  rolling  velocity  (p)  is  zero  by 
definition.  This  restraint  corresponds  to  0  =  ^  sin  6;  however,  for  small 
amplitudes,  the  two  axis  systems  are  comparable.  Although  errors  re¬ 
sulting  from  the  two -axis  systems  approach  are  believed  to  be  negligibly 
small,  in  general,  their  use  is  an  undesirable  characteristic  of  the  fitting 
technique  as  is  the  velocity  ratio  restraint. 

The  numerical  integration  fitting  technique  (NIF  technique)  developed 
for  axisymmetric  bodies  (see  Refs.  3  and  4)  is  significant  in  that  it  permits 
removing  some  of  the  undesirable  features  of  the  conventional  aeroballistic 
range  technique.  In  the  NIF  technique  for  axisymmetric  bodies,  the  differ¬ 
ential  motion  equations  are  referenced  to  a  fixed-plane  axis  system,  the 
velocity  ratio  restraint  is  removed,  and  nonlinear  aerodynamic  force  and 
moment  characteristics  can  be  handled.  An  additional  feature  of  this  tech¬ 
nique  is  that  the  motion  histories  of  more  than  one  flight  can  be  fit  simul¬ 
taneously.  This  multi-fit  feature  extends  the  capability  of  defining  non- 
linearities  permitting  the  use  of  flights  of  the  same  configuration  but 
which  experienced  different  mean  amplitude  levels.  However,  a  funda¬ 
mental  and  undesirable  feature  still  is  involved  with  the  NIF  technique  for 
axisymmetric  bodies  --  appreciable  aerodynamic  or  inertia  asymmetries 
can  not  be  handled. 

The  present  study  has  been  directed  toward  extending  the  capability 
of  the  NIF  technique  such  that  free-flight  motion  histories  of  bodies  with 
appreciable  asymmetries  can  be  analyzed  adequately.  Definition  of  the 
motion  equations  for  asymmetric  bodies  follows  the  same  approach  em¬ 
ployed  with  axisymmetric  bodies  in  which  the  basic  relations  of  Newton 
are  used: 


Here,  L  and  H  are  the  linear  and  angiilar  momentum  of  the  body  and  F 
and  M  are  the  corresponding  force  and  moment. 
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The  problem  of  defining  the  motion  equations  is  simplified  for 
asymmetric  bodies  by  using  a  fixed -body  axis  system  (rolling  axis  system), 
rather  than  the  fixed-plane  axis  system,  in  conjunction  y^ith  an  earth-fixed 
axis  system  (see  Fig.  lb).  The  time  rate  of  changes  of  the  linear  and 
angular  momentum  for  the  general  constant  mass  body  can  be  written 
(see,  for  example.  Ref.  6): 


Lx  = 

(u  +  qw  -  rv)in 

(1) 

S  = 

(v  +  ru  —  pw)ni 

(2) 

K  = 

(w  +  pv  -  qu)m 

(3) 

'xP  -  Ixy^  -  Ixz^  +  - 

V  -  P  q*xz  -  q^Iyz  +  +  ^P^xy 

(4) 

lyq  -  ‘xyP  -  'yz^  +  - 

K)  -  ^‘llxy  -  r^>xz  +  P^Uz  +  P  'Jlyz 

(5) 

•z"  -  IxzP  -  >yzq  +  P'i^Iy  “ 

‘x>  -  p^*xy  -  P^Iyz  *  I^Ixy  +  V  Ixz 

(6) 

a.  (X, ,  y, ,  I,  [  Frx^plane  axis  systam.  ^  ^ 

(xo.  yo»  *o)  Earth-fixad  axis  systam  ggrth-fixad  axis  systam 


Figure  1.  Axis  systems. 


As  all  current  and  proposed  reentry  type  missiles  have  at  least 
one  plane  of  symmetry,  the  angular  momentum  equations  are  simplified 
by  restricting  the  equations  to  bodies  for  which  the  (x-z)  plane  is  a  plane 
of  symmetry. 
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Hence, 

-H*  =  IxP  -  ^  ^  -  V  ■ 

Hy  =  lyq  -  pKl^  -  I,)  +  -  r^) 

Hz  =  Iz^  +  pq<iy  -  i;  +  Ixztq*-  -  p) 

Note  that  only  one  cross  product  of  inertia  component  (Ixz)  is  involved 
for  bodies  for  which  the  (x-z)  plane  is  a  plane  of  symmetry. 

Although  the  inertia  terms  associated  with  an  asymmetric  body  are 
easily  defined,  the  components  of  the  aerodynamic  forces  and  moments 
along  the  three  axes  of  the  body -fixed  axis  system  can  be  very  difficult 
to  define. 

In  Sketch  1  below,  the  angle  (a)  defines  the  orientation  of  the  plane 
of  the  total  yaw  angle  (6)  relative  to  the  body-fixed  axis  system.*  In 
contrast  to  the  usual  derivations  of  motion  equations  for  asymmetric 
bodies  (see,  for  example.  Ref.  6),  it  is  apparent  that  the  general  yawing 
motion  of  an  asymmetric  body  can  involve  in-plane  and  out -of -plane 
forces  and  moments  that  are  functions  of  a.  Of  the  forces  and  moments 


6  Cn 


(7) 

(8) 
(9) 


Sketch  1 
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that  are  dependent  on  the  roll  orientation  of  the  total  yaw  plane,  it  is 
believed  that  the  following  ones  are  of  primary  concern: 

1.  The  in-plane  and  out-of-plane  static  forces  and 
their  corresponding  static  moments  and 

2.  The  increment  involved  in  the  static  roll  moment. 

The  additional  aerodynamic  terms  corresponding  to  the  above  forces  and 
moments  are  the  ones  that  are  examined  here.  The  problem  of  defining 
the  aerodynamic  terms  for  an  asymmetric  body  can  be  better  appreciated 
by  first  examining  the  forces  on  the  axisymmetric  body  in  Sketch  2 
below; 


Sketch  2 


Again,  6  is  the  total  yaw  angle  and  a  is  the  orientation  angle  of  the  total 
yaw  plane  relative  to  the  y  and  z  axes.  From  aerodynamic  symmetry 
considerations,  there  will  be  no  out-of-plane  force  involved,  and  the 
normal  force  coefficient  in  the  5 -direction  can  be  written  for  a  body 
having  a  cubic  force  variation. 


sin  5  + 


where 


sin  S  = 
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and 


Hence, 


0 


^^N^S-DIR 


and  the  component  in  the  of- direction  is 


C 


N 


z 


cos  a 


•where 


cos  er 


2 


V 


Hence 


Replacing  and  Cj^gg  with  and  Cj^g^,  respectively,  in  the  last 
expression  follows  from  the  force  derivative  of  an  axisymmetric  body 
being  independent  of  a.  The  force  coefficient  in  the  ^-direction  can  be 
written  similarly.  The  point  of  interest  that  can  be  observed  in  the 
above  expressions  is  that  the  normal  force,  in  the  6 -direction,  for  an 
axisymmetric  body  having  a  linear  or  nonlinear  force  coefficient  is  a 
fimction  of  only  6.  Similarly,  in  the  case  o^*  the  linear  portion  of  the 
force  coefficient,  a  component,  say  in  the  a-direction,  is  dependent 
only  on  the  a-motion.  However,  the  or-component  of  the  nonlinear  portion 
of  the  force  coefficient  is  dependent  on  both  the  or-motion  and  the  /3-motion 
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experienced  by  the  body.  Similar  comments  obviously  hold  for  compo¬ 
nents  in  the  0 -direction.  It  is  important  to  recognize  that  the  total  normal 
force  for  an  axisymmetric  body  being  independent  of  a  permits  one  to 
easily  define  expressions  for  the  required  components  in  the  a-  and  i3- 
directions. 


In  the  case  of  an  asymmetric  body  having  a  linear  force  variation 
(see  Sketch  1),  Cn  and  Cno  are  the  corresponding  force  derivatives  if 
all  of  the  yawing  motion  is  Ather  in  the  a-  or  3-planes,  respectively. 

As  Cno  can  be  different  from  Cn  *  it  follows  that  normal  force 

derivafive  in  the  plane  of  the  tot^yaw  angle  (6),  will  necessarily  vary 
with  a,  the  orientation  angle  of  the  plane  of  the  total  yaw  angle.  Further, 
the  out-of-plane  force  derivative  [(CNg)o.  P.  I  ^  function  of 

cr.  It  follows  that  the  variations  of  Cjj  °  and  (CNg)Q  p  with  ct  are  re¬ 
quired;  however,  there  is  a  very  limited  amount  of  experimental  data 
available  to  aid  in  defining  precisely  the  variation  of  either  Cjq.  or 
(CNg)o  p  •  Hence,  approximate  variations  for  both  Cjijg  and  iCn5)o  p 
witn  0-  have  been  deduced.  These  are  given  below  and  are  believed  to" 
be  quite  adequate  for  bodies  for  which  the  force  derivatives  don't  change 
appreciably  between  positive  and  negative  yaw  angles. 

It  is  important  to  note  that,  in  the  defined  variation  for  C^gitwo 
basic  parts  are  involved.  One  part  is  dependent  on  the  roll  orienta¬ 
tion  (cr);  whereas  the  other  part  is  independent  of  ct.  The  approximate 
variation  for  used  is 

^Ns  =  ^ 


Here,  (CN5)jnean  portion  of  CN5  that  is  independent  of  <7.  Ob¬ 

serve  that  the  direction  of  ACn^  is  indicated  in  Sketch  3.  Also  shown 
in  Sketch  3  is  ACn^,  which  is  defined  here  to  be  the  difference  between 
Cn«  and 


Note  that  (CN5)jjjgan  corresponds  to  Cn^  in  the  ACj^jS  expression,  and 
that  ACn^  (dependent  on  ct)  is  a  function  of  AC^ig  and  can  be  written 


ACjyjg  =  ACjj  |sincr|  =  AC 


"/S  ”  /J7T 
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3 


y 


Sketch  3 


This  variation  for  ACjj^  follows  from  the  boundary  conditions  (see  Sketch 
3),  =  0  at  a  =  0  and  ACjij^  =  at  ct  =  90®.  The  expression  for 

can  now  be  written: 

%  -  y,2  (10) 

Further,  the  assumed  variation  of  (CNg)o,  p.  (see  Sketch  3)  with  a  is 


C 


o.p. 


Kjlsin  2a\  =  2Ki 


/~2  2 
Vw  +  V 


(11) 


This  variation  follows  from  the  boundary  conditions  [(Cn^)q  p  =0]  for 
a  =  0  and  90  deg  and  has  a  maximum  value  at  a  =  45  deg.  Note  that  Kj  = 
P.  at  CT  =  45  deg. 

The  resiilting  normal  force  expressions  derived  here  for  asym¬ 
metric  bodies  are  listed  below  and  are  consistent  with  the  following 
restraints.  The  portion  of  the  normal -force  coefficient  (linear  or  non¬ 
linear)  that  is  independent  of  a  has  been  handled  as  described  in  the 
previous  section  concerning  axisymmetric  bodies.  Portions  of  the 
normal -force  derivatives  that  are  dependent  on  a  have  been  assumed  to 
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vary  linearly  with  6 .  Hence  these  components  have  been  handled  con¬ 
sistent  with  the  previous  section  concerning  in-plane  and  out-of-plane 
force  derivatives  that  are  functions  of  <7.  Note  that  any  part  of  the  force 
derivative  that  is  nonlinear  and  dependent  on  a  has  been  assumed  negligibly 
small.  This  assumption  is  believed  to  be  quite  valid.  It  is  felt  that  these 
approximate  expressions  for  the  force  variations  with  a  are  reasonable, 
and  attempts  to  refine  them  further  are  not  justified  until  more  experi¬ 
mental  data  for  asymmetric  bodies  become  available.  Thus, 


The  corresponding  moment  expressions  including  terms  resulting 
from  the  above  in-plane  and  out-of-plane  forces  can  be  written: 


and 


C„  =  L 


+  2K2 


'-49' •  3-1 


C™  =  (■^)  - 

m  m  »  ^ 


+  2K. 


•  '-4;/  •  #] 


(14) 


(15) 


Finally,  for  missile -type  vehicles  which  usually  experience  small  roll 
rates,  the  u-dependent  increment  to  the  roll  moment  can  have,  in  some 
cases,  the  most  significant  o-dependent  effect  on  the  motion  of  the  body. 
The  assumed  variation  of  this  increment  with  ct  is 


ACg  = 


sin 


sin  5 
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This  variation  follows  from  the  boimdary  conditions  (see  Sketch  1), 

ACjg  =  0  for  <T  =  0  and  90  deg.  Hence,  for  the  present  elliptic  bodies, 

N  =  2,  and  the  rolling -moment  coefficient  as  used  in  this  report  can 
be  written: 

c,  -  (16) 

By  using  the  above  listed  expressions  for  inertia  terms  and  static 
aerodynamic  force  and  moments  and  also  accoimting  for  the  gravitational 
force  and  damping  moments,  the  basic  motion  equations  as  revised  for 
asymmetric  bodies  can  be  written,  viz. , 

P  =  Iz- V 


+  +  IzqAdCgi  2 

^xz 

(17) 

q=  [+pr(l2-I,()  + 

y 

(18) 

i  =  i-qrlx,(ly-lx-lz^  +  P‘i*L  -"Ix^Ix-VJ 

(19) 

+  Cj  +  IjjqAd  C^l - - 

I  V 


+  Cji  cos  0  cos  ^ 

(sin  d  cos  ^  sin  <^-  sin  ^  cos  0) 
+  Cjq  (sin  6  cos  \j/  cos  <f>  +  sin  ff/  sin  ^)] 


V  =  -  ^  V^[C  cos  0  sin 

®  Zm  * 

+  Cu  (sin  6  sin  \f/  sin  <f>  +  cos  jp  cos  <^) 

y 

+  Cj^  (sin  6  sin  if/  cos  <j)  —  cos  ^  sin  0)] 
z 

w  =  —  V^[—  C  sin  d  +  C«  X  sin  ^  cos  6 

”  2m  * 

-I-  Cjg  COS  ^  cos  6]  +  g 
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where. 


((»■  •  ?) 


^  V(w^  +  v^)  ^  ^^35 


Note  that  the  above  angular  motion  equations  define  the  motion  rela 
tive  to  the  fixed-body  axes,  whereas  the  translational  motion  equations 
define  that  motion  relative  to  earth-fixed  axes. 


The  programming  required  to  incorporate  the  motion  equations  for 
asymmetric  bodies  into  the  NIF  technique  was  accomplished  on  contract 
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■with  Armament  Systems  Department,  General  Electric  Company, 
Burlington,  Vermont.  Programming  aspects  of  the  NIF  technique  are 
discussed  in  detail  by  Whyte  and  Hathaway  of  the  General  Electric  Company 
in  Ref.  4.  The  fitting  programs  are  of  an  iterative  type  in  which  adjust¬ 
ments  to  the  equation  parameters  being  determined  are  obtained  using  the 
method  of  least  squares  and  differential  corrections.  These  programs 
have  considerable  flexibility  in  that  the  individual  aerodynamic  parameters 
listed  in  the  above  equations  may  be  fitted  or  held  constant.  Further,  the 
programs  include  additional  aerodynamic  terms,  not  listed,  which  permit 
handling  motion  histories  for  bodies  having  higher-order  nonlinearities  or 
for  the  more  simple  bodies  having  aerodynamic  and  inertia  symmetry. 

Some  of  the  roll -dependent  aerodynamic  terms  listed  in  the  above 
equations  reflect  modifications  to  terms  initially  included  in  the  present 
programs.  These  modifications  are  currently  being  incorporated  into 
the  programs. 


3.0  APPARATUS 


3.1  RANGE 


Range  G  consists  of  a  10-ft-diam,  1000-ft-long  tank  that  is  con¬ 
tained  within  an  underground  enclosure  (see  Fig.  2).  It  is  a  variable - 
density  aerodynamic  range  and  contains  53  dual-plane  shadowgraph 
stations.  Forty-three  of  the  stations  are  positioned  at  nominal  20 -ft 
intervals,  yielding  an  840-ft  instrumented  length.  The  spacing  of  the 
other  stations  is  such  that  a  station  is  located  approximately  10  ft 
downrange  of  stations  5  through  10,  12,  13,  15,  and  16.  The  angular 
orientation  and  position  of  most  test  configurations  can  be  determined 
to  within  approximately  ±0.  25  deg  and  ±0.002  ft,  respectively,  at  each 
station.  A  chronograph  system  measures  intervals  of  flight  time  to 
within  ±2  x  lO"*^  sec.  The  range  pumping  system  provides  range  pres¬ 
sures  from  1.  7  atm  down  to  about  20  juHg.  The  nominal  operating 
temperature  of  the  range  is  15°F.  The  launcher  used  in  this  investi¬ 
gation  is  a  two-stage,  light-gas  gun  having  a  2.  5-in.  -diam  launch  tube. 

3.2  CONES  AND  TEST  CONDITIONS 

Four  elliptic  cones  of  the  same  configuration  were  flight  tested. 

A  configuration  sketch  is  shown  in  Fig.  3.  The  cones  had  a  minor-to- 
major  axis  ratio  of  0.  75  and  were,  constructed  of  aluminum  with  a 
Fansteel®  60  nose  tip.  Two  pins  (0.025-  and  0. 04 -in. -diam)  were 
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inserted  into  the  base  of  each  cone.  These  pins  protruded  from  the 
base  0.  1  in.  and  provided  a  means  for  obtaining  the  roll  orientation  of 
the  cone  as  a  function  of  flight  time. 


Nosetip  Radii  = 

la  025 

(Nominal  4-percent  Bluntness) 


All  Dimensions  in  Inches 


Figure  3.  Elliptic  cone. 
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A  typical  four -component  sabot  was  utilized  in  launching  the  cones. 
The  sabots  were  constructed  of  Lexan®,  and  all  cones  were  launched  in 
an  uncanted  orientation  relative  to  the  sabot.  The  initial  angular  distur¬ 
bances  to  the  cones  were  those  arising  from  muzzle  effects  and  from  the 
cone -sabot  separation  process. 

The  four  cones  were  tested  at  a ’Mach  number  of  about  7.  5  and  a 
Reynolds  number,  based  on  cone  length  and  free -stream  conditions  of 
about  2.  7  X  10^.  The  physical  characteristics  and  test  conditions  of  the 
individual  cones  are  listed  in  Table  1 . 


Table  1.  Test  Conditions  and  Physical  Characteristics  of  the 
Elliptic  Cones  Investigated 


Shot 

Rej 

nip 

Iv  X  10*. 

lyXlO^ 
in. -lb-sec* 

Iz  X  10^. 

a 

No. 

M 

X  10-® 

gm 

In. -lb-sec* 

in.-Ib-sec^ 

a- plane 

1 

7.3 

2.6 

105.38 

1.08 

9.  11 

8.  79 

0.026 

0.033 

2 

7.4 

2.  6 

106.07 

1.09 

9.  23 

8.  89 

0.028 

0.  036 

3 

7.7 

2.8 

99.95 

1.09 

9.  32 

8.  97 

0.033 

0.039 

4 

7.6 

2.  7 

99.30 

1.07 

9. 29 

8.95 

0. 035 

0.047 

■Nose  Bluntness  Ratio 


4.0  UNCERTAINTIES  IN  AERODYNAMIC  PARAMETERS 


It  is  significant  to  note  that  the  experimental  errors  of  concern  in 
ballistic  range  aerodynamic  measurements  are,  in  general,  of  a  random 
nature.  Therefore,  it  is  felt  that  the  spread  in  the  measurements  of  an 
aerodynamic  parameter  provide  the  best  estimate  of  the  total  uncertainty 
in  that  parameter. 

Although  experience  in  reducing  data  using  the  current  fitting  pro¬ 
grams  are  limited,  uncertainties  in  the  reduced  parameters  should  be 
comparable  to  the  estimates  listed  below;  these  estimates  are  based  on 
spreads  in  measurements  observed  in  the  present  test  and  in  previous 
Range  G  test  programs. 

Coefficient 

Estimated  Total  1.  5 

Uncertainty  (20)*,  % 


2 


'N 


C£ 

0.  3 


'Ul, 


10 


-Standard  deviation  of  measurements 
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5.0  DISCUSSION 


Representative  plots  of  the  yawing  and  rolling  motions  of  the  elliptic 
cone  shots  are  shown  in  Fig.  4.  The  reversal  in  the  rolling  motion  (0)  is 
of  particular  interest  because  such  a  reversal  is  normally  not  experienced 
by  an  axisymmetric  body.  Satisfactory  fits  were  obtained  for  the  motions 
of  the  four  cone  shots  although  residuals  of  fit  were,  in  general,  about  25 
percent  higher  than  residuals  normally  observed  in  motion  fits  of  axisym¬ 
metric  bodies.  Accounting  for  the  roll  moment  contribution  that  is  ct- 
dependent  was  shown  to  be  very  important  as  the  inclusion  of  the 
term  was  necessary  to  permit  satisfactory  fits  of  the  motions  of  these 
shots.  However,  only  linear  aerodynamic  force  andyawing  moments  asym¬ 
metries  were  used  (Cn^  and  and  CmQ,  It  is  felt  that  any 

nonlinearities  associated  with  the  small  nose  bluntness  ratios  of  the  pre¬ 
sent  tests  would  be  negligibly  small;  however,  the  o -dependent  force  and 
moment  terms  discussed  previously  and  currently  being  incorporated  into 
the  fitting  programs  would  be  expected  to  permit  some  improvement  in 
the  motion  fits.  It  is  important  to  note  that  attempts  to  fit  the  measured 


ShotNg  1  Note:  Different  Symbols  Used  Only  for 

Unflagged  Symbols:  Curve  Rt  Points  Indicating  Different  Portions  of 

Flagged  Symbols:  Corresponding  Experimental  Points  the  Motion  Pattern 


deg 

a.  Yawing  motion 

Figure  4.  Representative  motion  patterns  for  the  elliptic  cones. 
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b.  Rolling  motion 
Figure  4.  Concluded. 


motion  histories  of  the  elliptic  cones  using  the  previously  developed 
fitting  programs  designed  for  axisymmetric  bodies  failed  completely. 

The  aerodynamic  data  obtained  for  the  elliptic  cones  are  presented 
in  Figs.  5  through  10.  In  Fig.  5,  the  components  of  the  normal  force 
coefficient  are  shown  as  functions  of  the  amplitude  experienced  in  the 
pitch  and  yaw  planes.  The  force  measurements  in  each  plane  are  con¬ 
sistent  for  the  four  shots.  For  the  cone  roll  orientation  indicated  in  Fig. 
5  (major  axis  of  the  elliptic  cross  section  in  the  cr-plane),  the  measure¬ 
ments  indicate  that  the  force  derivative  in  the  of-plane  is  approximately 
70  percent  of  that  measured  in  the  jS-plane,  which  is  indicative  of  an 
appreciable  configuration  asymmetry.  Motion  fitting  attempts  involving 
nonlinear  force  terms  did  not  produce  any  detectable  improvements  in 
the  quality  of  fits.  A  comparison  of  the  present  normal  force  data,  as 
a  function  of  ellipticity,  with  Newtonian  predictions  from  Ref.  7  is  shown 
in  Fig.  6  and  indicates  good  agreement  between  the  predicted  and  experi¬ 
mental  values. 
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(Without  Trim  Terms) 


Shot  No. 


Indication  of 
Maximum  Amplitude  in 
that  Plane  (Typical)- 


versus  a 
— —  Cify  versus  p 


e.  Cn, 


As  Viewed 
from  Uprange 
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The  moment  data  are  presented  in  Fig.  7  and  indicate  similar  effects 
of  the  configuration  asymmetry  as  observed  in  the  force  data.  Apparent 
deviations  in  the  measurements  for  shots  3  and  4,  and  for  shots  1  and  2, 
in  the  or-plane,  are  believed  associated  with  differences  in  nose  bluntness 
ratio,  (0)  (see  Table  1).  In  Fig.  8,  the  center  of  pressure  (cp)  data  are 
shown  as  a  function  of  the  nose  bluntness  ratio.  Though  the  measurements 
do  not  indicate  a  detectable  difference  in  cp  between  the  two  planes,  they 
do  suggest  a  bluntness  effect.  A  similar  effect  of  ^  on  cp  has  been  ob¬ 
served  previously  for  axisymmetric  bodies  (see  for  example.  Ref.  8). 

Note  that  the  observed  differences  in  nose  bluntness  resulted  from  con¬ 
struction  deviations  and  were  not  intentional. 


a,  3,  dee 


Figure  7.  Moment  coefficients. 
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a  68 
o 
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u 

S.  67 


u  66 


■Sym 

O  □  O  O  Different  Cones 
p-Plane 

Flageed  Symbols:  a-Plane 


z 


65  I _ I _ I _ I _ I _ I _ I 

0.02  0.03  0.04  0.03 


Nose  Bluntness  Ratio,  ^ 

Figure  8.  Variation  of  center  of  pressure  with  bluntness. 
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The  mean  damping  derivative  (Cmq)  is  shown  in  Fig.  9.  No  meas¬ 
urable  asymmetry  in  the  damping  moments  between  the  pitch  and  yaw 
planes  was  detected.  The  scatter  in  the  measurements  is  believed  related 
to  the  small  amplitude  levels  of  the  present  shots. 


-8 


O 


o  o 
o 


I _ I _ I _ 1 _ I _ 1 _ L 

6  7  8  9 

M 

Figure  9.  Damping  derivative. 


10 


The  mean  drag  coefficient  (C^)  is  shown  in  Fig.  10  as  a  function  of 
the  mean  amplitude  experienced  in  the  flights.  Attempts  to  detect  differ¬ 
ences  in  the  drag  coefficient  for  the  a-  and  ;3-planes  were  unsuccessful. 


6 

Figure  10.  Mean  drag  coefficient. 
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The  experimental  data  presented  are  significant  in  that  they  are 
believed  to  be  the  first  free -flight  data  obtained  for  such  bodies  and  they 
indicate  that  the  present  motion  analysis  procedure  provides  a  useful 
extension  to  the  free -flight  testing  capability.  It  should  be  noted  that  one 
can  expect  asymmetries  in  the  drag  and  damping  moment  always  to 
be  difficult  to  detect  at  small  amplitudes.  This  follows  because  the  effects 
of  these  two  parameters  on  motion  histories  are  small  when  small  ampli¬ 
tudes  are  involved. 

A  logical  extension  to  the  present  effort  would  be  further  revisions 
in  the  aerodynamic  expressions  to  account  for  magnitude  changes  in  a 
derivative,  such  as  Cmcy,  as  the  angular  motion  changes  sign.  Such  ex¬ 
tensions  are  obviously  desirable  as  they  would  permit  handling  asym¬ 
metries  of  greater  complexity. 

The  elliptic  configuration  of  the  present  test  has  two  planes  of 
symmetry,  and  hence  no  cross  product  of  inertia  terms  were  involved. 
However,  a  more  general  asymmetric  configuration  has  only  one  plane 
of  symmetry,  and  consistent  with  the  equations  of  this  report,  Ixz  can 
have  a  non-zero  value.  As  Ixz  be  difficult  to  define  experimentally,  it 
represents  an  area  of  concern  in  testing  general  asymmetric  configurations. 
It  was  initially  believed  that  there  was  a  possibility  that  1^^^  could  be 
treated  as  one  of  the  unknowns  in  the  motion  fitting  procedure  and  it 
covild  be  determined  similarly  to  the  unknown  aerodynamic  parameters. 
Hence,  in  the  present  programs,  the  motion  equations  were  programmed 
to  permit  an  option  of  treating  Ixz  either  a  known  or  an  unknown.  A 
limited  study  using  generated  motions  (with  and  without  simulated  ex¬ 
perimental  errors  in  the  model  orientation  angles  and  displacements) 
has  indicated  that  Ixz  determined  satisfactorily  in  the  fitting 

procedure.  To  date  simulated  vehicles  with  an  inertia  ratio  (Ixz/^y^ 
up  to  0.  06  have  been  examined,  and  the  aerodynamic  parameters  were 
evaluated  quite  well  when  Ixz  treated  as  an  unknown  in  the  motion 
fits.  It  is  important  to  note  that  for  the  higher  end  of  the  (Ixz/Iy)  ratio 
band  examined,  the  aerodynamic  parameters  could  not  be  determined 
if  Ixz  "was  ignored  in  fitting  the  motion,  that  is,  assumed  to  be  zero. 
Coirfirmation  that  Ixz  general,  be  handled  as  above  requires 

examination  of  a  wide  range  of  motion  patterns. 


6.0  CONCLUDING  REMARKS 

A  numerical-integration,  mot  ion- fitting  procedure  previously 
developed  for  axisymmetric  bodies  has  been  extented  to  vehicles  having 
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appreciable  aerodynamic  and  inertia  asymmetries.  The  primary  prob¬ 
lem  in  the  extension  of  the  procedure  is  related  to  the  required  revisions 
to  the  expressions  defining  the  aerodynamic  forces  and  moments  which, 
for  asymmetric  bodies,  can  be  very  dependent  on  the  roll  orientation  of 
the  body  relative  to  the  plane  of  the  total  yaw  angle. 

\ 

Experimental  free -flight  data  are  presented  for  a  cone  of  elliptic 
cross  section  that  has  appreciable  aerodynamic  asymmetries.  These 
data  indicate  the  potential  usefulness  of  the  developed  motion  an^ysis 
procedure. 
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NOMENCLATURE 


A  Reference  area  (cone  base  area) 

Cj3  Mean  drag  coefficient  for  flight 

Cjg  Rolling-moment  coefficient 

Cjg^  (T- dependent  rolling  moment  derivative 

Cjjj  Moment  coefficient  (a-plane) 

Cjnq  Mean  damping  derivative  for  flight,  ^  ^  (d^^v)  ^ 

Linear  and  nonlinear  moment  coefficients  (a-plane) 
Cmgj  ^msfi  Linear  and  nonlinear  moment  coefficients  (5 -plane) 
CivT  Normal-force  coefficient  (a-plane) 

Cj^  Normal-force  coefficient  (/3-plane) 

y 

0^3^  Linear  and  nonlinear  force  coefficients  (a-plane) 


^Nj3»  ^N3j3 


Cnj3»  CnsjS 


Linear  and  nonlinear  force  coefficients  (jS-plane) 
Linear  and  nonlinear  force  coefficients  (6 -plane) 
Moment  coefficient  (/3-plane) 

Linear  and  nonlinear  moment  coefficients  (3-pLane) 
Center  of  pressure  (percent  of  cone  length) 

Moment  reference  length  (cone  base  diameter) 
Angular  momentum  of  the  cone 


^xz»  ^xy> 


Moments  of  inertia 


Products  of  inertia 


Linear  momentum  of  the  cone 
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M 

m 

Rejg 

P,  q,r 

u,  v,w 

x,y,  z 
Xq,  yo*  Zq 
V 
a 
0 
6 
or 

0,  0 

SUBSCRIPTS 

x,y,  z 
0 

SUPERSCRIPT 


Mach  number 
Mass  of  the  cone 

Reynolds  number  based  on  free-stream  conditions 
and  cone  length 

Angular  velocity  components  of  the  fixed-body 
axis  system 

Linear  velocity  components  of  the  fixed-body 
axis  system 

Fixed-body  axis  system  (see  Fig.  1) 

Earth-fixed  axis  system  (see  Fig.  1) 

Velocity  of  cone 

Components  of  total  angle  of  attads  in  the  pitch  plane 
Components  of  total  angle  of  attack  in  the  yaw  plane 
Total  angle  of  attack 

Roll  orientation  angle  of  the  plane  of  the  total  yaw 
angle  relative  to  the  fixed -body  axis  system 
(see  Sketch  1) 

Euler  angles,  ^  is  also  nose  bluntness  ratio  (nose 
radius /base  radius) 


Indicate  directions  along  fixed -body  axes 
Indicates  coefficient  at  zero  yaw  angle 


Derivative  with  respect  to  time 
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